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I. INTRODUCTION 



The central problem in nuclear physics remains to understand the origin and nature of 
the nuclear force. In spite of the belief that we have attained the fundamental theory for the 
strong interactions- Quantum chromodynamics (QCD), this theory still eludes a satisfactory 
and complete description. The basic problem of QCD is that its natural and fundamental 
degrees of freedom, quarks and gluons, are not the observable baryon and meson states 
of the strong interaction. Thus bridging the missing link between the fundamental and 
observable degrees of freedom stands as one of the stark challenges of nuclear/elementary 
particle physics today. Although we do have an ab initio approach to solve this problem, 
that is lattice QCD, this endeavor is still miles away from achieving such a goal. This 
naturally motivates us to resort to non-perturbative QCD-based approaches of which this 
study is one. 

In the present paper, we address this lingering missing link by deriving a chiral meson- 
nucleon Lagrangian from a microscopic model of quarks and diquarks using path-integral 
methods. Chiral symmetry and its spontaneous breaking have consistently proven to be key 
concepts in understanding meson and baryon structure and many features of the nuclear 
force H], 0, ||, [|] . The gist of this paper is as following: We start from a QCD-based effective 
field theory to describe quark dynamics where the gluons have been integrated out. This is 
the SU(2)x, x SU(2)r Nambu-Jona-Lasinio (NJL) model that accommodates most of QCD 
symmetries [||, |J. Guided by general principles, we then assume that the nucleon can be 
described as quark-diquark correlations and introduce diquarks |7|, § as elementary fields in 
the problem. This assumption hinges upon the dynamical fact that two quarks can combine 
to form a color anti-triplet leading with the third quark to the formation of a color-singlet 
bound state, a baryon. Moreover, this assertion is vindicated by a mounting experimental 
evidence that diquarks play a dynamical role in hadrons |9|, [H| [IT], [12], [13], [14], [l^, [U| [17[| . 

We verify that only two kinds of diquarks are relevant for nucleons: the scalar isoscalar 
and the axial-vector isovector. By introducing composite meson and nucleon fields through 
the method of path-integral hadronization and then using a loop and derivative expansion 
of the resulting quark/diquark determinants, we arrive at an effective chiral meson-nucleon 
Lagrangian. The path-integral hadronization used here consists of two steps: bosonization to 
produce mesons as quark-antiquark correlations and what we label as "fermionization" which 
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generates baryons as quark-diquark correlations. In our model, mass, coupling constants, 
electromagnetic radii, anomalous magnetic moments, and form factors of the composite 
nucleon are calculated in terms of at most two free parameters. 

In this fashion, our treatment parallels, in the sense of calculating nucleon physi- 
cal observables, the approach of using the Faddeev equation |TS|, [Tj| for three quark 
, |2T], [22], P5| , or the approach of using static quark exchange |23J], the Salpeter 



states 



equation [^, |27|] or the fully relativistic Bethe-Salpeter equation |2^, ^9], |30[ for a quark- 
diquark system. Nonetheless, our formalism yields, in addition to nucleon observables, a 
Lagrangian of the quantum hadrodynamics (QHD) type [[31 , 



that describes the rich 



meson-nucleon interactions in a fully covariant and chirally symmetric formalism. 

While this program is applied to the case of deriving an effective Lagrangian for nucleons 
and mesons, it is certainly of general nature and can possibly be applied alternatively to 
yield prolifically other baryons and their interactions such as the A particle. Moreover, the 
idea of using path-integral techniques to transform a Lagrangian from its fundamental to 
its composite degrees of freedom is a powerful concept in physics of immense impact and 
utility. As a matter of fact, the authors of Ref. |33|] have recently invoked such path-integral 
techniques in their study of high-temperature superconductivity. They succeeded in doing 
so by converting a model of strongly-correlated electrons into an effective U(l) gauge field 
theory in terms of composite fields. 

The use of path-integral hadronization to derive a meson-baryon Lagrangian has been 
introduced in Ref. |34|, |35|, |36[ and applied to baryons with heavy quarks |37|, |38fl . Based on 
these ideas, the authors of Ref. |39|] attempted to construct such an effective Lagrangian for 
the nucleon using only scalar diquarks. They derived correctly the structure of the meson- 
nucleon Lagrangian, proved the Goldberger-Treiman relation and attempted to evaluate the 
axial-vector coupling constant qa as an application of their formalism. Their analysis and 
numerics for contain, however, few problems as well as an uncertainty due to the lack 
of a proper gauge-invariant regularization scheme. In the present paper, we extend their 
work by deriving the structure of the corresponding Lagrangian using both axial-vector and 
scalar diquarks, employ a gauge-invariant regularization scheme throughout our analysis, 
and verify the Ward-Takahashi identity and the Goldberger-Treiman relation. Furthermore, 
we present a full numerical study of various nucleon observables for the case of scalar diquarks 
drawing special attention to the role of an intrinsic diquark form factor. We concentrate our 
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analysis first on the scalar-diquark case for simplicity and due to the predicted dominance 
of this type of diquark in the nucleon ]35], fH], [LL], |42|] . Thus after more than ten years since 



the introduction of the idea of path-integral hadronization, this formalism is finally used to 
speak itself in calculating nucleon structure and its observables. 

The paper has been organized as follows. In Sec. II a microscopic model for quarks, 
diquarks, and their interactions is developed and meson and nucleon fields are introduced 



as auxiliary fields in the problem. The hadronization method is then invoked in Sec. [ITJ to 
rewrite the microscopic Lagrangian in terms of composite meson and nucleon fields. Next, 
a loop and derivative expansion is employed to calculate several terms in the Lagrangian 
including the nucleon self-energy and electromagnetic vertex. The issue of regularization is 
also examined and the Ward-Takahashi identity and the Goldberger-Treiman relation are 
verified. In Sec. a full numerical study for the nucleon is presented. Finally a summary 
and conclusions are provided in Sec. [V| as well as a discussion of some of the challenges and 
opportunities that remain. 

II. A MICROSCOPIC MODEL OF QUARKS, DIQUARKS AND THEIR INTER- 
ACTIONS 

A. Nambu-Jona-Lasinio Model 

In our model we treat the quarks using the NJL model which is a successful effective 
field theory where quarks interact through a four-point local fermion-fermion coupling. The 
highlights of the model are its incorporation of all global symmetries of QCD as well as its 
prediction of many features of QCD such as dynamical chiral symmetry breaking and its 
restoration H, 01 H3L H4|. Moreover, this model has been motivated, if also not derived, using 



lattice QCD g§, continuum QCD g3} fPj, g§, gU|, and Yang-Mills theories 0. The 



locality assumption of the model has been justified for low energy QCD [|50f| , and inspired 
by strong-coupling lattice quantum electrodynamics (QED) |9|, |5Tfl . The main problem of 
the NJL model continues to be the absence of confinement. Therefore, the success of the 
model rests on observables that are insensitive to the details of confinement. It is noteworthy 
here that there exists various attempts to include the effects of confinement within the NJL 



model [43, 52 1 . 
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We start from an NJL Lagrangian satisfying SU(2)l x SU(2)^ chiral symmetry: 

G 

Cnjl = q(ifi - m )q + — [(qq) 2 + {qi^fqf] . (1) 

Here q is the current quark field, r are the isospin (flavor) Pauli matrices, G is the NJL 
coupling constant, and mo is the current quark mass which explicitly breaks chiral symmetry. 
The color and flavor indices are suppressed in this expression and assumed to be so for the 
rest of the paper unless explicitly shown. Starting from this Lagrangian, we construct the 
corresponding vacuum partition function as 

Z=M l J VqVq exp % J d 4 x [£ NJL ] , (2) 

where M\ is a normalization constant. 

B. Introduction of meson fields 

Composite scalar (a ~ qq) and pseudoscalar (7? ~ qij^rq) meson fields are introduced as 
auxiliary fields in the problem. This is done by multiplying the NJL partition function of 
Eq. (0) by the term (with J\f 2 being another normalization constant) 

1 



M 2 I VaVn exp i / d 4 x 



2G (^) 



(3) 



At this stage, no modifications have been made to the underlying dynamics of the Lagrangian 
as this multiplicative factor is merely an overall constant in the partition function. We 
impose the following transformation: 

a > a + Gqq , 

7r' ; — > 7r* + Gqirf^q , (4) 

in order to eliminate the quadratic terms (~ {qq) 2 ) of the NJL Lagrangian. Using transla- 
tional invariance of the integration measure T>aT>7T, this results in the expression: 

1 



W 3 J VaVnVqVq exp i J d 4 x 



q(ip-m -a - i^r ■ 7?) q - T^O 2 + ^) 



(5) 



The prescribed change in field variables is nothing but the Hubbard- Stratonovich trans- 
formation j53, 54]. We label the resulting Lagrangian as the "semi-bosonized" one since 



we have already introduced the boson (meson) fields but have not yet integrated over the 



quark ones. The current quark mass mo is then absorbed into the definition of the field a 
and the meson fields are further transformed according to the non-linear parameterization 

[<7,7r]-[<7',$] 



a + ry 5 f ■ vf = (m q + a') exp ( -— 75T ■ $ ) , 



(6) 



44 1 . Accordingly, the NJL 



where = 93 MeV is the pion decay constant and m q = (a)o is the constituent quark mass 
which is fixed through a gap equation in the meson sector 
Lagrangian is converted to 

1 



ift - (m q + a') exp ( --^f ■ $ 



£ N jl = 5£ sh - ^ (</ + m q f + q 



Here, the 8C S \> = (9 (mo) is the symmetry-breaking mass term given by 

i 



(7) 



5C 



sb 



m q + a 
16G 



mo tr 



exp 



-7 5 r • $ +h.c 



(8) 



where the trace is taken over flavor and Dirac indices. 



C. Diquarks 

In studying the Lorentz and flavor structure of the qq correlations, we find five possible 
types: the scalar qq, pseudo-scalar qi^^q, vector qj^q, axial-vector qj^j^q, and tensor qa^ u q 
diquarks. Here q = q T C~ x ^iT2 where T stands for transpose and = ^7270 is the inverse 
charge-conjugation matrix. Moreover, we identify two isospin structures for each of these five 
Lorentz qq formations. Explicitly, we have an isoscalar and isovector diquarks by inserting 1 
and r between q and q. It has to be noted here that the q spinor has the same transformation 
properties as q in the Lorentz and isospin groups. 

A question arises concerning how many of these ten diquarks are needed to form the 
nucleon. Using permutation symmetry and Fierz transformation, we have verified an earlier 



assertion |55| that only two diquark formations are independent for the nucleon if the nucleon 



field is to be written as a local operator of three quarks. This result is consistent with the fact 



that in constituent quark models [56], the nucleon wave function is constructed using only 
scalar and axial-vector diquarks. Hence, we introduce two diquarks as elementary complex 
fields: as an axial-vector isovector field with electric charge = {4/3, 1/3, -2/3} and D 
as a scalar isoscalar field with a charge = {1/3}. 



D. Quark-diquark interaction terms 



Quark-diquark interaction terms are introduced to form the nucleon as a relativistic 
bound state of quarks and diquarks. We consider such an interaction in a local form. 
This is essentially the static approximation of solving the three-body equations for baryons 



within the NJL model |22| , 42]. It is more convenient here, in terms of forming a chirally 
invariant quark-diquark couplings, to work with the chirally rotated "constituent" quark 
field x defined by 1 



1 ,5 T 



X = exp f -—7 5 - • $J q . (9) 

The range of possible symmetry preserving interaction terms is limited 2 . This provides a 
highly welcomed dynamical constraint in our treatment. Discarding for a moment interaction 
terms describing a possible scalar-axial- vector mixing (see Eq. (|i~3|) below), we may choose 
the following term for the quark-scalar-diquark interaction: 

C qD ~ G ( X D^) (D X ) , (10) 

while we may select 

£ q m ~ G (xYl'r- Dt) (A, ■ f tVx) , (H) 

for the quark-axial-vector-diquark coupling. Our choice for the full interaction term is 
dictated by the need to produce the nucleon as a linear combination of axial-vector and 
scalar diquarks according to 

B ~ G (sine D^-t^J 5 X + cosO D X ) ■ (12) 

In the above expressions G is the quark-diquark coupling constant with mass dimension 
[G] = m _1 , while 9 is a mixing angle for the two diquark contributions. These are, as we 
shall discuss below (see Sub-Sec. [IV A| ), the only free parameters in our model. 



1 The chiral rotation of Eq. (||) induces anomalous terms as the phase of the integral measure. In this 
paper, we are concerned only with the non-anomalous processes. 

2 If we work with the current quark field g, that is in the linear representation of chiral symmetry, a 
vector-isoscalar diquark is necessary as a chiral partner of the axial- vector-isovector one. 
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E. Microscopic Lagrangian 



Electromagnetic interactions are introduced in our model through the canonical method 
of covariant derivatives in the quark and diquark Lagrangians. We proceed to form the 
microscopic Lagrangian by batching the diquark contributions, the quark-diquark interaction 
terms including mixing, and the semi-bosonized NJL Lagrangian of Eq. ([?]) after the field 
transformations of Eq. (ffi), to obtain the following Lagrangian as our input model: 



£ = xS-'x - ^' + m ^ + + D^A-'D + D ] M K~l D v + 

G (smO xYr 'f f - &l + cos^x^ 1 ) (sin D v ■ r 7 Vx + cos9 D X 



(13) 



where 



s- 1 

M 



+ M , 

Y\ ■ V; + -A; + a' + >fQ q A 



EM 



(14a) 
(14b) 



Here S 1-1 is the modified inverse propagator that includes the free inverse propagator Sq = 
(iP — m q ) for the quark field and the interaction matrix Ai. Note that Q q = diag (2/3, —1/3) 
is the quark charge. The M. matrix contains all interaction vertices of the quark with meson 
and electromagnetic fields. The quark interacts with the pion through vector V£ and axial- 
vector A* functions of the pion field. These arise from the derivative term ifi following 
the transformation of Eq. (pi). Precisely, these functions are defined through the Cartan 



decomposition (£ = -^-): 

exp ( ~1*T ■ £ j exp ( % -^f ■ £ 



l -l 5 r-AlH) + |r-V;(0- 



(15) 



The M. matrix also encompass pion-photon and weak gauge boson vertices which are not 
shown for brevity. 

In the above Lagrangian, we use the modified scalar diquark inverse propagator 



A' 1 = A^ + iQ s Af A {P -P) 



(16) 



and the modified axial-vector diquark inverse propagator 

A" 1 = Ao - 1 + iQ A \{A™K - ATX) - 9,uAT(P - P) 



(17) 
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where we have omitted 0(Q 2 S ) and 0{Q\) terms for brevity. Each of these expressions 
includes the free inverse propagator as the kinetic and mass term, and the electromagnetic 
interaction vertex. Explicitly, the free inverse propagator for the scalar diquark is 

A" 1 = -(«9 2 + Mf) , (18) 

while the one for the axial-vector diquark is 

K-l = 1® [g^(d 2 + M\) - d,8 v ] , (19) 

where 1 is the unit matrix in the three-dimensional isospace. Here Qs = 1/3 and Qa = 
diag (4/3, 1/3, — 2/3) are the charges of the scalar and axial-vector diquarks respectively. 
The modified axial-vector diquark propagator encompasses (not shown) weak interaction 
terms. 



F. Introduction of nucleon fields 



There is still one missing component in our microscopic model: collective nucleon fields. 
Therefore, a nucleon field B is introduced as an auxiliary one by multiplying the partition 
function of the Lagrangian of Eq. ( |T3D by the term 



M A J VBVB exp i J d 4 



x 



-Ibb 

G 



and 



where A4 is a normalization constant. In a similar fashion to Eq. 
transform the field configuration according to: 

B — ► B + G (sin9 D u -rYl 5 X + cos9 D X 
B — ► B + G (sin 9 xYl 5 r ■ t>\ + cos9 

As a result, the quark-diquark interaction term in Eq. flT3] ) is rewritten as 

^J- BB - B (sin9 D u -rYl 5 X + cos9 D X 
G ^ 



(20) 
), we 



(21) 



A/4 / VBVB exp \ i d 4 x 



- (sin 9 X7 M 7 5 r- + cos9 xD^ B X . (22) 

This procedure completes the introduction of composite meson and nucleon fields into the 
problem and wrap up the construction of the microscopic model. 



III. DERIVATION OF A MESON-NUCLEON LAGRANGIAN 



A. Hadronization of the microscopic model 

At this point we have a Lagrangian that involves only quarks and diquarks as dynamical 
fields with kinetic and mass terms while the meson and nucleon fields are merely auxiliary 
ones. Additionally, the quark and diquark fields appear in bilinear forms appropriate for 
integration as a consequence of eliminating the interaction terms. Thus, we rearrange the 
expressions involving the quark fields into the form \ S~ X X ~~ VX~ XV) an d use the fermion 
path- integral identity (Det stands for determinant): 



/ V X Vx exp i J (xS-\ - f K - xn) = Det (5"') exp (-i / f, S „) 



(23) 



to integrate over the quark fields. In doing so, we would have accomplished the path-integral 
bosonization that delivers to mesons their full dynamical character || |44|. We still need to 
integrate over the axial- vector and scalar diquark fields in order to achieve a meson-nucleon 
Lagrangian. Thus, we cast the terms that involve the diquark fields into the form tp'tCtp 
where (p = ^-D M , D j , and use the boson path-integral identity: 

J V^Vpexpi J (ipJC^) = [Det(/C)] _1 . (24) 

This final integration procedure is what we label as fermionization as it produces fermions 
from boson-fermion correlations. The quark-diquark dynamics has been absorbed by the 
composite meson and nucleon fields. We have at last fully "hadronized" the quark and 
diquark Lagrangian. The microscopic model of quarks and diquarks has been converted into 
a "macroscopic" model of mesons and nucleons possessing the same (approximate) chiral 
symmetry as the original microscopic fields. Notice that the quarks and diquarks do now 
appear only as virtual particles in loops and are described by corresponding propagators 
and interaction vertices. 

Next, we use the relation Det(M) = exptrln(M), to rewrite the determinants as La- 
grangian terms of meson and nucleon fields. Thereupon, we arrive at a compact chiral 
meson-nucleon Lagrangian given by 

C cS = 5C sh - -L(a' + m 9 ) 2 - ztrlnS- 1 - I BB + i tr ln(l - □) + 
2(jt G 

i tr ln(l - A EM Int) + i tr ln(l - A EM Int) . (25) 
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Here the trace is over color, flavor, and Lorentz indices while the "EM Int" label stands 
for the electromagnetic interaction terms of each of the diquarks as given in Eq. ( pp and 
Eq. fllTD- Furthermore, 




(26a) 



where 

S = cos 2 6 B A S B , (26c) 

= sin 9 cos 9 B A S r j j v j 5 B , (26d) 

(jr 2 )Mi = sin q cos q Q A pk >^ 7p7 5 r k S B . (26e) 

The configuration and color indices have been omitted for simplicity. The effective hadron 
Lagrangian of Eq. (|25D contains plenty of rich physics. It encompasses, through the loop 
and derivative expansion in the tr In terms, kinetic and mass terms for nucleons and mesons 
together with a multitude of possible interaction terms of mesons, nucleons, and electroweak 
gauge bosons. It comprises terms describing the various electroweak interactions of mesons 
and nucleons such as meson photoproduction (the Kroll-Ruderman terms ||57||), and in ad- 
dition, it includes terms for meson-meson and nucleon-nucleon scattering. Nonetheless, the 
most desired part of the Lagrangian is the prized meson-nucleon interaction and nucleon- 
nucleon vertices which delineate the nuclear force. 



B. Self-energy diagram and kinetic terms 

The physics in the Lagrangian becomes manifest in terms of loop and derivative ex- 
pansions of the resulting quark-diquark determinants. We concentrate here on the nucleon 
sector which is contained in the terms % tr ln(l — □) — A BB. Nicely, the expansion 

tr ln(l -□) = -tr (□ + — + — + •.. J , (27) 
is an expansion in the number of nucleon fields. Moreover, since 

S = (1 + SoMy^o = S -S MS + So M So M S + • • ■ , (28) 
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B(p) 



o 



B(p) 



p-k 



FIG. 1: The Feynman diagram for the nucleon self-energy which generates the nucleon kinetic 
and mass terms and produces the mass equation that determines the nucleon mass. 

each term in the logarithmic series leads to an expansion in terms of the number of interaction 
vertices : 0, 1, 2, 3, • ■ ■ . 

We take the first term in the two expansions (Eq. (p7|) and fl2"8|)) which is up to integrations 
(see Fig. 0) 



B(x) 



E(x,y) + -^8{x-y) 



B{y) 



(29) 



where S(x,y) is the nucleon self-energy: 

E(x,y) = —N c % 7 " 7 B r< iA*^(ar,y) iS (x,y) Tj y»y 5 sm 2 8 
-N c iA (x, y) iS (x, y) cos 2 6» . 



(30) 



Here A^ c = 3 is the number of colors (resulting from the trace over color), and summation 
over repeated indices is to be understood. The Fourier transform S(p) of the self-energy is 
decomposed according to 



E(p) = S s (p 2 ) + ^E,(p 2 



(31) 



The nucleon mass Mb is then given by the vanishing of the inverse nucleon propagator: 



- + E S (M|) + M B Z v (M 2 B ) = 0. 

Gr 



(32) 



This condition generates dynamically the nucleon mass (one of the predictions of the model) 
in terms of the theory parameters, and it is similar in structure to the mass equations that 



determine the masses in the meson sector || 44J]. Near the mass shell, the inverse nucleon 



propagator takes then the form: 



E(p) + ^ 



(J) — Mb ) Z 



-i 



(33) 
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Evidently, the nucleon has now acquired the desired status as a dynamical degree of freedom 
in the problem. Here Z is the wave-function renormahzation constant (see Sec. |A 1| ) which 
prompts us to renormahze the nucleon field according to B = y/Z B vcn . 

C. Regularization of divergent integrals and Ward identity 

Now we are in a place to discuss regularization. The self-energy and the various Feynman 
diagrams in the problem involve the evaluation of divergent integrals. Consequently, we are 
confronted with the question of how to regularize these integrals. This issue emerged as a 
decisive one in our analysis as we have attempted several regularization schemes. We started 
by adopting the method of four-momentum sharp cut-off, but found it unsatisfactory as it 
yielded a violation of the Ward-Takahashi identity. Guided by "experience" , one can remove 
by fiat the terms that violate gauge invariance, and thus conforms to this identity |39 



However, a more rigorous and solid method is certainly desirable. Accordingly, we sought to 
regularize the integrals using both the three-momentum sharp cut-off and the Pauli-Villars 
methods. The former is motivated by dispersion theory ||3| and leads, as we verified, to 
compliance with the Ward identity. Yet, we found that the most suitable regularization 
scheme, in terms of rigor and convenience, is the Pauli-Villars technique which we have 
established in this work as the standard method for regularizing all divergent integrals. 
This method consists of introducing a fictitious propagator with some mass M to cancel 
the divergent contribution in the integral at large momentum values (see the Appendix |A|). 
As a matter of principle, the Pauli-Villars mass which appears in the nucleon sector can be 
different from the NJL cut-off A arising in the meson sector. Nonetheless, to minimize the 
number of free parameters, we elected to equate them, M = A. It is noteworthy here that all 
observables (see Sec. [IV| below) were found to be insensitive to the value of the Pauli-Villars 
mass upholding the futility of using it as a free parameter. 

In the process of testing gauge invariance (the Ward-Takahashi identity), we have to 
determine the electromagnetic vertex of the nucleon. This implies evaluating two kinds 
of diagrams depicted in Fig. || where the nucleon can couple to the electromagnetic field 
through either the quark or the diquark propagators. For the Ward-Takahashi identity 
to be satisfied, the wave-function renormahzation Z must be equal to Z\ where Z\ is the 
electromagnetic vertex renormahzation constant at q 2 = 0. Here, q = p' — p (incoming 
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FIG. 2: The Feynman diagrams for the electromagnetic coupling which generate the nucleon 
electromagnetic vertex. 




FIG. 3: The Feynman diagram for the axial- vector coupling which generates the quark contribu- 
tion to the nucleon weak-interaction vertex. 

photon) is the momentum transfer. This condition is indeed satisfied for both the Pauli- 
Villars and the three-momentum cut-off methods. By calculating these diagrams at an 
arbitrary value of momentum transfer, we derive the nucleon form factors from which we 
can extract the electromagnetic radii and anomalous magnetic moments. 



D. Nucleon axial- vector vertex and the Goldberger-Treiman relation 

We are in a position to calculate the weak-interaction axial-vector vertex (Fig. |3|) to 
determine the axial- vector form factors from which we can extract the axial- vector coupling 
constant g&. The Qa is defined as the coefficient of the 7^7 5 term of the axial- vector vertex 
at vanishing momentum transfer. This vertex leads naturally to the Goldberger-Treiman 
relation as can be seen by noticing that 

= i d,$ + 0($ 3 ) . (34) 
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Thus we obtain the following term for the pion-nucleon coupling: 

g A B Yl 5 ^ ■ $1 B — > j^- B -f-? T - -d^B. (35) 

But this term has to be identified with the pseudovector form of the Yukawa pion-nucleon 
coupling |3"2|| : 9 *£ N B ^^ b f- B, which prompts us to conclude that 



g<s>NN = -et9a ■ (36) 

This is nothing but the Goldberger-Treiman relation at the composite hadron level. Note 
that this relation appears intact with no effort in our treatment as opposed to large violations 
of up to 30% in the Bethe-Salpeter equation approach p9fl . 

Finally, one must mention that in addition to the weak-gauge-boson coupling through the 
quark line (Fig. |3[), there is, only for the case of the axial- vector diquark, a weak coupling 
through the diquark line. The quantitative contribution of this diquark to the Goldberger- 
Treiman relation remains to be investigated in a future work. 

IV. NUMERICAL STUDY 

Having thus far derived the structure of the problem, we proceed to generate numerical 
results for our model using the simpler case of only scalar diquarks, thereby admitting the 
possibility of an intrinsic diquark form factor. Including only scalar diquarks is not out 
of place. Indeed, recent studies using scalar diquarks have reported good results for most 



of the nucleon observables p4| HI- Moreover, there are indications of scalar diquark 
dominance in the nucleon [5], [IT]]. As a matter of fact, a very recent calculation using 
the Faddeev equation for three quark states has concluded that axial-vector correlations, 
while still important for magnetic properties, contribute at most no more than 10% to the 
structure of the nucleon ^ . 



A. Free parameters and basic quantities 

Tab. | provides the basic quantities in our model. As for free parameters, we have first 
the NJL coupling constant G and the NJL cut-off A which are fixed to yield the constituent 



quark mass (quark condensate) through the NJL gap equation in the meson sector [EL H3, [44 
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TABLE I: Basic quantities in our model: The basic quantities in the microscopic model are the 



constituent quark mass 


(mq), the scalar diquark mass 


(Ms), the cut-off A, and the quark-diquark 


coupling constant G. 






m q 


M s 


A G 


.390 GeV 


.600 GeV 


.630 GeV 271.0 GeV" 1 



In this fashion, the G constant decouples completely from the nucleon sector, the sector of 
our interest. The two diquark masses are also determined, in a consistent manner, using 
the NJL model and the Bethe-Salpeter equation in the diquark channels^, OT. In this 
context, the diquark masses are simply poles, just as mesons, but in the quark-quark T 
matrix. This leaves us with only one new free parameter in our model: the quark-diquark 
coupling constant G. As can be discerned, this model is well-constrained and yields a 
powerful predictive strength. It has to be remarked here that in principle there is another 
free parameter in the model : 9 as the mixing angle for the two diquark contributions. But 
this angle has no effect in the present analysis as we consider only scalar diquarks (6 = 0). 
Moreover, one must mention that while the quark and diquark masses and the cut-off A are 
in principle fixed through the NJL model, small variations in their values are permissible as 
they still lead to consistent results within the NJL model. This adds a margin of freedom 
to these masses. 



B. Nucleon static properties 

Tab. H displays our predictions for some of the static properties of the nucleon. Experi- 
mental values are taken from Ref. |58| , |59|j . A mass of 0.94 GeV is obtained for the nucleon 
through the mass equation Eq. (|3^). By fixing the nucleon mass at this value, we would have 
eliminated the G coupling constant from the problem and reached a theory with no more 
free parameters. As for the binding energy of the nucleon, it is estimated (m q = 0.390 GeV 
and M s = 0.600 GeV) as AE h - m = m q + M s — M B = 50 MeV, suggesting that the nucleon 
is a loosely bound state of a quark and a diquark. This prediction is consistent with other 
approaches of the NJL model or using the Faddeev equation |^2| . 



In the same table, we show the magnetic moments of the proton and the neutron. Our 
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TABLE II: Model predictions: Some of the nucleon static properties as predicted in the present 
calculation using the intrinsic diquark form factor (IDFF) or without it. Experimental values are 
taken from Ref. 58, 59]. 





f' P 


fin 


9A 


<r 2 > P E 


< r 2 > E 


< r 2 > p 


< r 2 










(fm 2 ) 


(fm 2 ) 


(fm 2 ) 


(fm 2 ) 


Theory with IDFF 


1.57 


-.75 


.87 


.77 


-.11 


.82 


.84 


Theory without IDFF 


1.57 


-.75 


.87 


.68 


-.19 


.82 


.85 


Experiment 


2.79 


-1.91 


1.26 


.74 


-.12 


.74 


.77 



treatment predicts a number that is two-third of the experimental value for the proton 
and about one-half of that for the neutron. This is not a surprising result considering 
that we have not included the axial-vector diquark in the present calculation. Constituent 
quark models and other more sophisticated approaches predict precisely that the axial- vector 
diquark inclusion should add the missing one-third strength to the proton and the missing 
one-half one to the neutron |24], p6fl . So in this context, this result is exactly what one should 
have expected from our current analysis. For the same reason, the predicted value for the 
axial- vector coupling g& of 0.87 is significantly less than the experimental one of 1.26. While 
the scalar diquark cannot couple to the weak interaction, the axial-vector one indeed does 
couple to the weak gauge bosons adding strength to the interaction. The magnetic moments 
and the axial-vector coupling constant display a rather small sensitivity to the coupling 
constant G. In fact the calculated values for the magnetic moments are not that different 



from the predictions based on simple additive models |2J] suggesting a rather independence 
from the details of the dynamics or the nucleon size. 

Speaking of the nucleon size, it is nicely well-produced by our model: the electric and 
magnetic radii for the proton and the neutron are close to the experimental measurements. 
The negative charge radius of the neutron has been suggested as an indication of a scalar 
diquark clustering in the nucleon |60[| , and our treatment manifests this conjecture in a 
dynamical model. These numbers point to a physical picture of a "heavy" diquark at the 
center with a quark rotating around it. By comparing the radii as calculated with and 
without the intrinsic diquark form factor, we find that the extended size of the diquark 
contributes a positive value of about 0.10 fm 2 for each of the proton and neutron electric 
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radii. Note that the scalar diquark has a positive charge and thus the contribution is positive 
adding about 0.10 fm 2 to the proton radius and reducing the absolute value of the neutron 
one by the same amount. We conclude here that the size of the diquark contributes about 
10% of the proton radius and 40% of the neutron one. This confirms an earlier calculation 
using a static quark-exchange approximation [|24|| . As anticipated, the intrinsic form factor 
has virtually no effect on the magnetic radii as the scalar diquark, as we shall see below, 
has a negligible contribution to the magnetic form factors. Moreover, we have found the 
electromagnetic radii to be very sensitive to the binding energy as well as, although implicitly 
through the mass equation (Eq.( ^)), the coupling constant G. This indicates, as easily 
expected, the importance of the details of the dynamics for the nucleon size. We should 
remark here that the quantities calculated here do not contain the pion contribution which 
becomes significant for some physical quantities. It is estimated that for m n = 138 MeV, 
typical values for the pion corrections are of order of 30% 3 . 



C. Nucleon electric and magnetic form factors 

Next we calculate the nucleon form factors. Fig. || displays the proton electric form factor 
in comparison with experimental data taken from Ref. ||61|| . The figure shows the quark 
contribution, the diquark contribution in addition to the full form factor (the sum of the 
two contributions), with no intrinsic diquark form factor. It is evident that our treatment 
reproduces beautifully the form factor at low values of momentum transfer (Q 2 = —q 2 
where q is the momentum transfer). The discrepancy at higher values of Q 2 begs for an 
understanding. The figure suggests an explanation: the diquark contribution is almost 
constant implying a rather localized diquark inside the proton. Although treated as an 
elementary field in our theory, the diquark is a composite object and does have a finite 
size. Our treatment needs to be adjusted to reflect the true nature of the diquark by 
incorporating an intrinsic diquark form factor. Fortunately, diquark form factors have been 
calculated recently and in fact in the framework of the N JL model |Sj . Thus we readily add 
the intrinsic diquark form factor to our treatment and produce the proton form factor shown 
in Fig. [|. Impressively, the calculations matches very well with the experimental data. 

3 In the exact chiral limit (m^ — > 0), the pion contribution to the isovector charge radius diverges. 
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FIG. 4: The proton electric form factor in comparison with experimental data. The figure shows 
the quark contribution (dotted line), the diquark contribution (dashed line) as well as the full form 
factor (dotted-dashed line) as the sum of the two contributions. Intrinsic structure of the scalar 



diquark is not included here. Experimental data from Ref. [61] are included 



The neutron electric form factor tells a similar story. In the left panel of Fig. |6|, just as in 
Fig. |], we display the neutron form factor with its quark and diquark contributions. Clearly, 
the quark contribution is negative in value (d-quark) and thus cancels much of the diquark 
contribution leading to a small form factor. There is once more a discrepancy compared 
to experimental data that is largely eliminated once we include the intrinsic form factor 
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0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 

Q 2 (GeV 2 ) 

FIG. 5: The proton electric form factor, with intrinsic diquark form factor (diquark f.f.), in 
comparison with experimental data. The figure shows the form factor with no intrinsic diquark 
form factor (dotted-dashed line), and with the intrinsic form factor (solid line). Experimental data 
from Ref. |3l| are included. 



as exhibited in the right panel of the same figure. It is noteworthy here that the neutron 
form factor is a potent test of any treatment as it is a delicate cancellation of two large 



contributions ]26], P7| . Saliently, the cancellation is naturally produced in our study. The 
experimental data are obtained from Ref. |4], |65|, |66|, |67| . 

In Fig. [7|, we present the proton magnetic form factor as calculated with or without the 
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FIG. 6: The neutron electric form factor in comparison with experimental data. The left panel 
shows the quark contribution (dotted line), the diquark contribution (dashed line) as well as the 
full form factor (dotted-dashed line) as the sum of the two contributions. The right panel displays 
the neutron form factor with no intrinsic diquark form factor (dotted-dashed line), and with the 
intrinsic form factor (solid line). Experimental data from Ref. (63, |64|, B^, |66|, B7| are included. 



intrinsic diquark form factor. The figure also contains the quark and diquark contributions. 
Unmistakably, the scalar diquark contribution is virtually vanishing due to the lack of an 
intrinsic spin. Nevertheless, there is a very small contribution due to a small orbital angular- 
momentum effect in the bound quark-diquark system. Since the diquark contribution is 
negligible, the inclusion of the intrinsic diquark form factor does not alter our prediction 
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and the form factor is determined to be almost purely from a quark origin. This suggests the 
need for the axial- vector diquark, which does have an intrinsic spin, to supplement the quark 
contribution and to provide the missing one-third strength compared to the experimental 
data [EH]. The figure also indicates a convergence of our calculation, mainly from a quark 
origin, and the experimental data at large values of Q 2 . Such result suggests that this form 
factor is almost purely from a quark origin in this regime. This is anticipated due to the 
finite size of the axial- vector diquark which probably can have a significant contribution but 
only for smaller values of Q 2 . 

The neutron magnetic form factor describes a similar narrative to that of the proton but 
here the missing strength (one-half) is larger as can be seen in Fig. ^| As stated earlier, 
these specific missing strengths are predicted due to the absence of the axial- vector diquark 
in the present analysis. The experimental data in Fig. ||] are obtained from Ref . |j4 , [68|. |69| . 



V. CONCLUSIONS 



In this paper we tackled the nucleon structure and the challenging problem of under- 
standing the origin and nature of the nuclear force by deriving a meson-nucleon Lagrangian 
using the path-integral method of hadronization. We started from a microscopic model of 
quarks and diquarks where the gluonic degrees of freedom have been integrated out. The nu- 
cleon was conceived as quark-diquark correlations and only two kinds of diquarks were found 
relevant for its structure. These are the scalar isoscalar and the axial-vector isovector di- 
quarks. Composite meson and nucleon fields were introduced by the methods of path-integral 
bosonization and fermionization to rewrite the problem in terms of the physical meson and 
nucleon degrees of freedom. This yielded an effective chiral meson-nucleon Lagrangian after 
using a loop and derivative expansions of the resulting quark/diquark determinants. The 
divergent loop diagrams were regularized using gauge-invariant regularization schemes and 
the Ward-Takahashi identity and the Goldberger-Treiman relation were verified. 

An extensive set of nucleon observables were calculated for the first time on the basis of 
the path-integral hadronization approach. Indeed, many of the nucleon physical properties 
such as mass, coupling constants, electromagnetic radii, anomalous magnetic moments, and 
form factors have been determined from a model of essentially one free parameter. By taking 
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FIG. 7: The proton magnetic form factor in comparison with experimental data. The figure shows 
the proton magnetic form factor as calculated with (solid line) or without (dotted-dashed line) the 
intrinsic diquark form factor. The figure also includes the quark (dotted line) and the diquark 
(dashed line) contributions to this form factor. Note that the three curves (apart form the diquark 



contribution) are very similar. Experimental data from Ref. [Sl| are included. 



into account the intrinsic diquark form factor, we established a remarkable agreement with 
the experimental data for the nucleon size and the electric form factors, while our calculations 
show missing strengths for the magnetic form factors and the axial- vector coupling constant. 
The discrepancy is likely due to the absence of the axial-vector diquark in the present 
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numerical study. 

This work is part of an ambitious program of using path-integral techniques and QCD- 
based effective field theories to study baryon structure and to derive a full-fledged nuclear 
force. The final goal of the program is a derivation of an effective field theory for the nuclear 
force, namely, quantum hydrodynamics (QHD) from quark dynamics. As for the future, we 
plan to attain a numerical study using both the scalar and axial-vector diquarks including 
their intrinsic form factors . This is a challenging and rather difficult task due to the axial- 
vector diquark intricate structure as a particle of one unit spin and isospin. Some of the 
ensuing complications are the axial-vector diquark direct coupling to the weak interaction 
and the electroweak scalar-axial-vector transitions. At a later stage, we plan to generalize 
our approach to chiral SU(3) symmetry to study the structure of the baryon octet. 



APPENDIX A: ANALYTICAL EXPRESSIONS USING THE PAULI-VILLARS 
REGULARIZATION METHOD 

We include here analytical expressions for some of the principal formulae in our treatment. 



1. Self-energy and wave- function renormalization 



The self-energy is depicted by the Feynman diagram of Fig. [I] and is given by the expres- 



sion: 



S(p) = N c 



d 4 k i {j> — jt + m q 



(2tt) 4 (p-k) 2 -m 2 



k 2 -M 2 k 2 - A 2 



N, 



1 



c (4tt) 2 



dx — x) + m q ] In 



A(A,g 2 = 0) 
LA(M s ,g 2 = 0)J 



(Al) 



Here p M is the momentum of the nucleon taken to be on the mass shell, N c = 3 is the number 
of colors, and 



A(M, q 2 ) = m 2 x + M 2 (l - x) - p 2 x(l - x) + q x 2 y(y - 1) . 



(A2) 



Note that we have used the Pauli Villars method for regularizing the divergent integral by 
incorporating the propagator k i^ A i of a fictitious scalar particle with mass A in the above 
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expression. 

The wave-function renormalization Z is obtained through the derivative 9 a ^ 



, lead- 



*m b 



ing to the expression: 
1 



Z 



(4vr) 2 









(1 — x) In 







A(A,g 2 = 0) 



2x(l - x) M B [M b (1 -x) + m q 
2. Electromagnetic interaction 



A(M s ,g2 = ) 
1 



1 



A(A,g2 = ) A(M s ,q 2 = 0) 



• (A3) 



The photon couples to both the quark and the diquark lines leading to two contributions 
to the nucleon electromagnetic vertex. 



a. Quark contribution 

This is the contribution represented by the left part of Fig. ^ and is given by: 

d 4 k 



^ q = ~N c Q q 



(2tt)< 



k 2 - Mj k 2 - A 2 



+ m q ) ^ i(^>-jt + m q ) 
(p' — k) 2 — m 2 (p — k) 2 — m 2 ' 



2 . io** q v 



2Mp 



(A4) 



where F^(q 2 ) and F|(g 2 ) are the quark contributions to the nucleon form factor and are 
given by: 

A(A,g) 
A(M s ,q) 



F^q 2 ) = ZN c Q q -^ J dxj^ dy { (1 - x) In 



2x(l - x)M B [M B (1 -x) + m q \ + gV(l - y) 



|_A(A, 9 ) A(M s ,q) 



F 2 V) = ZN c Q q — — j dx j ih, 



x 2x 2 M B [Mb(1 - x) + m q 



(A5) 



(A6) 



|_A(A, q) A(M s ,q)\ 

In the above expressions, p M (p ,fl ) is the momentum of the incoming (outgoing) nucleon, 
p^ is the momentum transfer, and Q q = diag (2/3, — 1/3) is the quark charge. 



P 



Note that we have included the wave-function renormalization constant Z in the above 
expressions. 
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b. Scalar diquark contribution 



This is the contribution represented by the right part of Fig. ^| and is given by: 

_A'' -—NO [ d * k f 1 i(jt + m q ) i 

7D ~ c ^ S J (2tt) 4 1 (p - k) 2 - M 2 k 2 -m\ {p'-k) 2 -M 2 



i(k + m a 



(p-k) 2 -A 2 k 2 -m 2 (p'-k) 2 -A 2 



(p" - p' - 2k") 



W) >f + F?{q 



Df 2^ ^ Qu 



2M, 



B 



where 



ZN C Q t 



1 



dx / dy 



(4vr) 2 , .0 
2x(l - x)M B [M B (1 - ar) + m g 



1 — x) In 
1 



A(A,g) 
LA(M 5 ,g) 



A(A,g) A(M s ,q) 



F 2 u (q 2 ) = ZN C Q, 



(4tt) 2 



i ,i 
dx / dy 

JO 



2x(l - x)M B [M B (1 - z) + m ff ] 
Here Qs = 1/3 is the diquark charge. 



A(A,q) A(M s ,q) 



(A7) 



(A8) 



(A9) 



c. Form factor and Ward-Takahashi identity 

The full nucleon electromagnetic form factor is the sum of the quark and diquark contri- 
butions. Note that at q 2 = 0, the sum of these two pieces conforms to the Ward-Takahashi 
identity as it yields Z (Q q + Qs) -§ = diag(l,0) = Qn with the correct normalization for 
the nucleon electric charge. 



3. Axial-vector coupling constant qa 

The axial-vector coupling constant Qa (with only scalar diquarks) is determined from the 
axial- vector vertex as represented by Fig ^. Note that since the scalar diquark cannot couple 
to the weak interaction, there is no contribution to Qa from a direct coupling to the diquark 
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line. The full axial-vector vertex is given by: 



M 1 = - N 

-^axial iV c 



d 4 k 



k 2 -M 2 k 2 - A 2 



i{f-jt + m q ) 5 i{j)-k' + m q ) 



(p' — k) 2 — m 2 (p — k) 2 — m 2 



(A10) 



The gA is defined as the coefficient of the 7 M 7 5 term of the axial-vector vertex at q 2 
This yields after evaluating this vertex 



g A = -Z N c 



(47T) 2 Jo J 



dx l dy < xln 



A(A,g 2 = 0) 
A(M 5 ,g 2 = 0) 



x [M B (1 — x) + m q 



A(A,g 2 = 0) A(M s ,g 2 = 0) 



• (AH) 
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FIG. 8: The neutron magnetic form factor in comparison with experimental data. The figure 
shows the neutron magnetic form factor as calculated with (solid line) or without (dotted-dashed 
line) the intrinsic diquark form factor. The figure also includes the quark (dotted line) and the 
diquark (dashed line) contributions to this form factor. Note that the three curves (apart form the 



diquark contribution) are very similar. Experimental data from Ref. [64, 68, f39], |7^, [71], [T^] are 
included. 
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